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ABSTRACT 



I find the common envelope (CE) energy formalism, the CE a-prescription, to be inadequate 
to predict the final orbital separation of the CE evolution in massive envelopes. I find that when 
the orbital separation decreases to ~ 10 times the final orbital separation predicted by the CE a- 
prescription, the companion has not enough mass in its vicinity to carry away its angular momentum. 
The core-secondary binary system must get rid of its angular momentum by interacting with mass 
further out. The binary system interacts gravitationally with a rapidly-rotating fiat envelope, in a 
situation that resembles planet-migration in protoplanetary disks. The envelope convection of the 
giant carries energy and angular momentum outward. The basic assumption of the CE a-prescription, 
that the binary system's gravitational energy goes to unbind the envelope, breaks down. Based on 
that, I claim that merger is a common outcome of the CE evolution of AGB and red super-giants 
stars with an envelope to secondary mass ratio of Mcnv/-^^2 ^ 5. I discuss some other puzzling 
observations that might be explained by the migration and merger processes. 



INTRODUCTION 



The common envelope (CE) process is in the heart of the formation of many close binary syste ms (e.g., 
Iben fc LiviolllQQi iTaam fc Sandauistlbood IPodsiadlowskilboOll IWebbinklliooj iTaam fc Ric keil boiol). During 



the CE phase the orbital se paration decreases due to gravitational drag and tidal interaction (e.g. Jlben fc Livio 
I993I : iRicker fc Taamll2012r) . The transfer of orbital energy and angular momentum to the envelope, as well as 
other possibl e energy sources , lead to the ejection o f the envelope. One of the major unsolved questions of the 
36 years old ( Paczvnski 1976t van den Heuvel 1976 ) CE process is the final orbital separation. 



In the commonly used energy formalism of the CE the grav itational energy released by the spir aling-in 
binary system Eq, is equated to the envelope binding energy (e.g.. IWebbinklll984l : iTauris fc Dewill200ll) . i?bind, 
with an efficiency a: aEo — i?bind- This is termed the CE a-prescription. It is now a c ommon practice to 
include the internal energy of the envelope in calculating the binding energy (e.g. iHan et al. |[l994; Maxted et al 



2002: Zorotovic et al 



2O10l:IXu fc Lilboid 



Davis et al 



2011), and some autho rs argue for observational support 



for that (e.g., Rebassa-Mansergas et al. 2012). Ivanova fc ChaichenetsI (2011) suggest that the enthalpy rather 
than the internal energy should be included in calculating the binding energy. 

The basic assumption of the CE a-prescription is that the binding energy is channelled to eject the envelope 
in a uniform manne r (for a thorough discussion of this prescription and other aspects of the CE evolution see 
Ivanova et aL 2012fl. Name l y, the re is no separation between envelope parts. This assumption was put into 
question bv iKashi fc Soken (|2011l ) in c ases where the final CE phase is a rapid process. While some studies 
do suggest rapid final evolution, (e.g., Rasio fc Livio 1996t Livio fc Sokei 1988 ), on the order of days to few - 



week s , others argue for a final phase that la sts for months (e.g., Sandauist et al. 1998t De Marco et al. 2003 



20091 : iPassv et al.ll2011t iRicker fc Taamll2012l ). In any case. iKashi fc Sokeij (j201l[ ) studied the rapid deposition 
of energy in the very inner regions of a massive envelope of an asymptotic giant br anch (AGB) sta r , and found 
that a substantial fraction of the ejected mass does not reach the escape velocity. iKashi fc Soken (|201lh used 
a self-similar solution and followed the blast wave propagation from the center of the AGB outwards. They 
showed that ^1 — 10 per cent of the ejected envelope remains bound to the remnant binary system. They 
further argued that the bound gas falls-back and forms a circumbinary disc around the post-CE binary system. 
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The interaction of the circumbinary disc with the binary system wiU reduce the orbital separation much more 
than expected in the CE a-prescription. The smaher orbital separation favors a merger at the end of the CE 
phase or a short time after, wh ile the core is still hot. A diffe rent mechanism for merger at the termination of 
the CE phase was suggested by Ivanova fc ChaichenetsI ( 2011 ). 



A suppo rt to the separat i on of t h e envelope to ejected an d bound segments comes from numerical simulations 
as well, e.g., Lombardi et al. ( 20061) . Sandauist et al. (1998) showed that some of the mass of the envelope can 
remain bound to one or both of the interacting stars. They simulated a 5 Mq AGB interacting with a 0.6 Mq 
companion, and found that the companion unbinds ~ 1.55 Mq (~ 23 per cent) of the AGB envelo pe. Th e y also 
obtained a differenti ally rotating thick disk or torus at intermediate stages of the CE evolution. Sokei ( 19921 ) 
(see also ISoken 12004) had analy t ically obtained a similar thick disk structure. In these cases a rapid merging 
is expected as well. iPassv et al.l (|201l[ ) present more extreme results in their simulations. They find that when 
the envelope is lifted away from the binary > 80 per cent of the envelope remains bound to the binary. In their 
simulations rotation of the AGB was not included, so the material that remains bound is probably overestimated. 
They conclude t hat in some case s parts of the AGB envelope remain bound or marginally bound to the remnant 
post- AGB core. iDe Marco et al.l (|201l[ ) suggest that envelope material which is still bound to the binary system 
at the end of the CE will fall back onto the system and will form a circumbinary disk. They suggest that such 
a disk might have some dynamical effects on the binary period. 



Ricker fc TaamI (|2012( ) performed 3D numerical simulation of a low mass binary composed of a 1.O5M0 red 
giant and a O.6M0 companion. Although the AGB e nvelope mass is much lower than the masses considered 
here, some results are relevant . Ricker fc Taam |201 j) find that during the rapid inspiral phase (the plunge-in 



phase; see Ivanova et al. 20121 ) only a fraction of ~ 25% of the energy released by the spiraling-in process goes 



toward eject i on of the envelope. The mass loss rate at the end of their simulation is ^ 2M0 yr ^. As mentioned, 
Passv et al.l (j201ll ) find that most of the envelope s tays bound at the end of t heir simulations (on the comparison 



between the simulations of Passv et al. 2011 and Ricker fc Taam 20121 see Ivanova et al. 2012). These results 



might imply that for a white dwarf (WD) companion of mass ~ O.6M0 inside an AGB star of mass > 4:Mq, the 
envelope can stay bound for few years after the orbital separation has shrunk to several solar radii. At the end 
of the simulation the envelope is concentrated around the equatorial plane of the binary system. This result of 
a rapidly rotating disk- like envelope structure is one of the motivation for the present paper. 

In the present paper I put into question some basic assumptions of the CE a-prescription. I do not dispute 
that the CE helps in ejecting the envelope. However, most of the gravitational energy is released near the 
final expected separation. There, I argue in section [51 some assumptions of the CE a-prescription break down. 
Despite that, the CE a-prescription might give a crude estimate of the final separation in cases that don't end 
up in merger. The final migration to merger is discussed in section [3J and a phenomenological toy model is 
proposed in section U) Summary and implications for some puzzling observations are discussed in section O 



2. INTERACTION DEEP IN THE COMMON ENVELOPE 

The density profile of AGB stellar envelopes can be approximated by p{r) oc r^^ where /3 ~ 2 — 2.4. For 
the purpose of this paper it is adequate to take for the envelope density 

where po is the density at a reference radius rg, M^nv is the envelope mass, and i?* is the stellar radius. The 
envelope mass inside radius r is given by 



The binding energy of the envelope is given by the expression 

£^bind = 0.5 / 47rr p[r)dr, (3) 
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where the factor 0.5 comes from the usage of the virial theorem and Afcoro is the core mass. Because the 
integration diverges at low radii for the simple model used here, the binding energy is calculated for the envelope 
residing above radius r„i ~ IRq. Performing the integration then gives the binding energy 



-Bbind(fm) 



— -1 



+ 



Mr. 



where the second equality defines F. For a typical values of 
Tm — l-R©, we find F ~ 9. 



cnv-^''^^corc -, 



(4) 



0.7 Mq, ~ 3OOi?0 and 



According to the CE a-prescription the final orbital separation would be determined by the equality 
0.5aGMcoroAf2/aQ — -Ebind, where M2 is the secondary mass and Oq, the final separation according to the 
a prescription. This gives 

F\-' 

(5) 



Oq M2 

— = O.Ola — 

R* 0.2Afciiv 



10 



Substituting in equation ^ gives the envelope mass inner to that radius 

-0.2 / \ -0.8 



MM ^ o_,3^o.8 
Mo 



A/2 



0.2Afp 



(6) 



For this scaling we find that A/er(2aa) = 0.22Af2 and Afer(3aa) = O.5M2, and that M^rir) = M2 at 
r = ISac,. This result implies the following for CE inside evolved red giant stars (AGB, RGB, etc.). As the 
orbital separation shrinks to ~ IOoq, the companion has not enough mass in its vicinity to carry away its angular 
momentum. Hence, the interaction with the gas in the vicinity of the companion cannot cause much further 
spiraling-in, and no gravitational energy is liberated t o unbind the envelop e. At this point the so called plunge-in 
phase, where the spiraling is dynamical and rapid (see Ivanova et al.l[2012t) . must end. The core-secondary binary 
system must get rid of its angular momentum by interacting with mass residing at r > a, where a is the orbital 
separation at the given moment. In other words, local interaction with the envelope gas at a typical distance 
of the Bondi-Hoyle-Lyttleton accretion radius and tidal interaction with the mass residing at r < a, cannot 
cause further spirallin g-in. The binary system interacts gravitationally with a rapidly-rotating flat envelope 
(jRicker fc Taamll2012l ). where flat envelope refers to a highly oblate envelope structure. This situation becomes 
more like the planet-migration process in protoplanetary disks rather than the hydrodynamical interaction of a 
secondary star in the classical prescription of a CE. At this stage the flat envelope dissipates the energy deposited 
by the spiraling-in binary system, and creates heat that can be transferred outward by convection. The basic 
assumption of the CE a-prescription, that the binary systems gravitational energy goes to unbind the envelope, 
breaks down. 

This discussion raises questions about the en tire CE a-prescrip t ion. P r oblems for the CE a-pre scription 
in explaining observations were noted before, e.g., Nelemans fc Tout ( 20051) . Nelemans fc Tout J 2005 ) find the 
CE-7 prescription to better fit observations. In the CE— 7 prescription (jNelemans et al.l 120001 ) the angular 
momentum released by the spiraling-in binary system is assumed to be carried away by mass loss. The CE— 7 
prescription also seems to break-down at a < lOoa, for the same reason the CE— a does: there is not enough 
mass in the vicinity of the secondary a nd inward to transfer the angular momentum to. In any case, the CE— 7 
prescription has other severe problems (jlvanova et al.ll2012l) . 



3. MIGRATION-INDUCED MERGER 

The conclusion from the above discussion is that when the orbital separation inside RGB and AGB stars 
decreases to a radius where M{a) ~ M2, which occurs at a ^ Cq, the spiraling-in process substantially slows 
down. Here as before, a is the binary separation and is the final separation predicted by the CE a-prescription. 
Although numeri cal simulations do not always find rapid envelope rotation or a very flattened envelope (e.g.. 



Passv et al.ll201lh . I assume here that after the slowing down of the spiraling-in process, such a state will be 



achieved. The binary system transfers angular momentum to the rapidly-rotating circumbinary flat envelope, and 
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consequently t he average radius of the e nvelope expands, and the binary separation decreases and the eccentricity 
increases (e.g.. lArtvmowicz et al.lll99l|). I assume that t he str ucture of the flattened e nvelope crudely resembles 



the structure of the disk studied in lArtvmowicz et al.l (|l99lh . In the calculations of lArtvmowicz et al.l (|l99ll) 
the disk extends from the nearest stable circumbinary orbit of ^ 2 .5a up to 6a. I assume that the binary orbit 



is circular (eo — 0). Scaling the results of lArtvmowicz et al.l (|l99ll ) by values appropriate for the present study 
gives the rate of change of the semi-major axis 



-4.5 X 10" 



(It) 



where 



is the binary orbital angular velocity, and 



2.7 



qd = 



M„.(6a) 



lAMp 







(7) 



(8) 



Mcorc + M2 

Here Mor(6a) is the envelope mass inner to r = 6a. The spiraling-in timescale at a = li?© with this scaling is 



(9) 



(10) 



The Keplerian (orbital) time decreases with decreasing separation, but then so does the mass in the flattened 
envelope. The time scale of a year is long enough to allow the convective envelope to transport the heat outward. 
Therefore, instead of expelling the entire envelope, the spiraling-in process increases the luminosity. The wind 
mass loss rate might increase substantially, but a large portion of the envelope will stay bound in systems where 
the initial envelope mass is Mcnv ^ several x M2. The binary system will merge. 

Equation Q, although crude for our case, none the less seems to be adequate for the present goals. In a 
very recent paper IShi et al.l (j2012i) performed 3D magnetohydrodynamic (MHD) simulations of a circumbinary 
disk surrounding an equal mass binary system. Similar to earlier studies, they find that strong torques by the 
binary clear a gap of radius ~ 2a. The binary system gains angular momentum by accretion from the disk, and 
loses it by the gravitational torque. Over all the orbital separation decreases. We note that the accretion rate 
onto the two stars in the present study is expected to be very low, or practically zero, due to the high pressure 
ar ound the c o re of the AGB star. To the contrary, mass is expected to flow outward. Although the torque found 
bv IShi et al.l (|2012r) is much larger than the one used in the ad— viscosity models (not to be confused with the 
CFi—a), because of the opposite effect of accretion they find the spiraling-in r ate to be only a factor of ^ 3 
higher than found here. The numerical coefficient in equation (O according to IShi et al.l (j2012l ) is —8 x lO^'^. 
The difference is of no significance for the goal of the present paper. 

For type II migration, where the planet is relatively massiv e and a gap is opene d in the disk, the migration 
rate is limited by viscous transport in the disk. The rate given bv lAlibert et al.l (j2005h . for example (their eq. 19), 
gives a similar value to that in equation ([7]) here when the ad— v iscosity coeffic i ent is taken to be aj, — 0.01 — 0.1. 
The same holds for the type II migration expression derived bv lNelson et al.l (|2000l ). for example. 



Over all, I conclude that the timescale for the final spiralling-in process for systems with Afonv ^ several x 
M2 is practically the viscous time scale in the flattened envelope. The major contribution to the viscosity is 
convection. This implies that on the same time scale that the binary system spirals in, the convection transport 
the released gravitational energy out. Instead of being used to expel the entire envelope, the released binary 
gravitational energy is radiated away. The increase in luminosity will substantially expand the envelope and 
increase mass loss fro m the surface. Howe ver, mass will stay in the vicinity of the binary system, and might 
lead to a merger (also lKashi &: Sokeiil201ll ). A different process, based on the r eexpansion of g a s above the cor e 
and inner to the secondary orbit, that might lead to merger was discussed bv llvanoval (j2011r ). Ilvanoval (|201ll ). 
as well as other papers dealing with post-CE merger, though, do use the energy formalism for the CE ejection 
(the CE a-prescription) . Here the fundamental assumptions of the CE a-prescription are criticized. 
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4. A PHENOMENOLOGICAL SCENARIO 

The termination of the CE phase with a core-secondary binary system interacting with material further out 
is hkely to be more violent than the migration of a planet (and more complicated) . Here I demonstrate a simple 
phenomenological approach as a demonstration of a plausible alternative to the CE a-prescription at the phase 
when the mass inward to the secondary is less than the mass of the secondary. At this preliminary stage this 
should be considered more as a toy model rather than an established process. 

I assume that bound mass is falling toward the binary system. No real disk is formed, but rather when 
the bound gas gets to a distance of rja, where 77 ~ 1, from the binary system it acquires energy and escapes. 
For simplicity consider the case where the core and secondary mass are about equal, Mcoro — If the 

secondary mass is much below that of the core, the interaction will not be violent, and migration, more like 
planet migration, takes place. Consider that the core-secondary system interacts with bound gas of mass Mb, 
part of which fell to the center, at a typical distance a from the center of mass. Each binary component is at 
a distance of a/2 from the center of mass. As the core and secondary come closer to each other, so does the 
surrounding bound gas. The binary system ejects the gas. Namely, the gas that is ejected carries a specific 
energy of deb — G(Afcore + M2)/(ria). The change in orbital separation is given by energy conservation 

, /^M^\ ^ g(Mco..e + M.)^^^^^ ^^^^ 
\ a J rja 

where Mb is the ejected mass from the binary vicinity (hence increases with time). Integration of the left 
hand side from the orbital separation ab at the time the circumbinary interaction phase starts to a final orbital 
separation af, and the right hand side from zero to Mb, yields, for A/corc = M2, 

af ~ abe-^''"/^'"^ = 06(0.018)*^'/"^'^* (12) 

where = Mcoro + M2. 

In this toy model the final orbital separation is very sensitive to the circumbinary bound mass present when 
the migration starts. I emphasize again that this toy model applies only when the secondary mass is not much 
smaller than the primary mass. For af, — IORq and rj ^ I, for example, the final orbital separation in this toy 
model is a/ = 3.7Rq, IARq, O.Si?©, and 0.2i?Q, for Mb/Mt — 0.25, 0.5, 0.75, and 1, respectively. The last two 
orbital separations might lead to a rapid merger. 



5. DISCUSSION AND SUMMARY 

In the previous sections I questioned the basic assumption of the common envelope (CE) energy formalism 
(the CE a-prescription) , and found it to be inadequate to predict the final orbital separation of the CE evolution 
in massive envelopes. The CE a-prescription might help as long as the envelope mass inward to the location 
of the secondary is larger than the secondary mass. This does not occur in the inner region when a massive 
envelope is involved, as evident from equation ©. Based on that, I claim that merger is a common outcome of 
the CE evolution of AGB and red super-giants stars with an envelope to secondary mass ratio of AIcnv/M2 > 5. 
Many binary systems do survive of course. Here Monv is the envelope mass at the beginning of the CE phase. 



Me rgers at the termination of the CE phase were discussed before (e.g., Ilvanova fc Podsiadlowskil 12003 



Ivanovalboill). as well as the slowi ng down of the spiraling-in process when the separation is small (see discussion 



Ivanova fc Podsiadlowskil 2003h . However, these authors and many others do use the CE energy formalism. 



In the energy formalism of the CE the binary gravitational energy that is released is equated to the envelope 
binding energy. Most of the binary gravitational energy is released very deep in the envelope. It was found here 
that where most of the gravitational energy is supposed to be released, there is not enough local envelope mass 
to absorb the energy released by the binary system (for AGB stars with Menv/M2 > 5). The energy must be 
absorbed by material residing further out in the flattened envelope, a process that proceeds on a viscous time 
scale, as with migration of planets inside disks. Viscosity is supplied by the convective envelope. The convective 
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speed in AGB stars is close to the sound speed, and the same convective motion can transport the energy out. 
A large fraction of the released energy will go to increase the luminosity and inflate the envelope, rather than 
to expel the inner parts of the envelope. 

Merger at the termination of the CE phase might explain some puzzling observations. I list here several 
of these, each of which deserve much de eper study. I do not list all possible outcom es of CE mergers, e.g., the 
formation of high-fi eld Magnetic WDs ( Tout et al. 2008t Garcia-Berro et al. 2012 ) and NS-BH mergers (e.g., 
Dominik et al.ll2012r ). but rather limit myself to those scenarios that are theoretically less developed. 



The Core-Degenerate (CD) scenario for Type la supernovae. ISoarks &: Stechen (|1974) developed a merger 
scenario of a WD with an AGB core, in which the result is a type II SN. In a later study Livio fc Riess ( 20031) 
found that the merger of a WD with an AGB core can lead to a SN la that occurs at the end of the CE phase 
or shortly after . Hyd roge n lines will app e ar, an d hence they attribute such a scenario to a rare type of SNe la. 
Kashi fc Sokeii (j201ll) and lllkov fc Soken (|2012l ) argue for a long time delay between the CE merger and the SN 
la explosion due to a long spin-down process, and termed it the core-degenerate (CD) scenario for SN la. They 
argue that this scena rio might account for most SNe la. An important conclusion of iLivio fc Riess ( 2003) an d 
Kashi fc Sokeii (|201ll ) is that for merger to occur the AGB star should be massive. As explained in lSokeii (|2011[ ). 
the CD scenario has some key differences from the double-degenerate scenario for SNe la, and hence should be 
referred to as a separate scenario. 

R Coronae Borealis stars with a massive expanding shell. The two leading models for the formation of R 
Coronae Borealis sta rs are a final-helium-shell flash of a post- AGB star and two WDs merger. In a recent paper 
Clavton et al.l (j201l[ ) find R Coronae Bo realis to have signatures that are ex pected from both a WD merger and 
from a final- helium-shell flash (see also iKameswara Rao fc LambertI 120111 ). The relatively high inferred mass 
of R Coronae Borealis and its high fluorine abundance support the merger model, while the massive, ~ 2M0, 
expanding shell supports the final-helium-shell flash model. I raise here the possibility that R Coronae Borealis 
and similar R CrB stars are formed from the merger of a WD with the core of an AGB star. 

The bright edge of the planetary nebula luminosity functio n fPNLF). The hright edge of the PNLF in old 
and young stellar populations is the same (for a recent paper see Ciardullol 2012 ) . Two scenarios were suggested 
to account for the puzzling finding that old stellar populations have as bright planetary nebulae (PNs) as young 
stellar populatio ns do. One is that the bright part of the PNLF in old populations comes from blue-straggles 
(|Ciardulld 120121 ). namely, merger of main sequence stars. The other one is that the brightes t PNs in old 
populations are systems in transition from a symbiotic nebula phase to the PN phase (Soker 20061) . Here I raise 
the following scenario to account for bright PNs in old stellar population. In a system of two low mass stars, 
the more massive lost its envelope to its companion while on the red giant branch (RGB). A helium WD of 
^ 0.2 — O.3M0 is formed. Later, when the secondary becomes an AGB star, the WD remnant of the primary 
star merges with the AGB core before the hydrogen-rich envelope is lost. The now more massive core can ionize 
the nebula to form a bright PN. This speculative scenario deserves more detailed study. 

The formation of single sdB stars. Extreme horizontal branch star (EHB; also termed sdB or sdO stars) 
are horizontal branch stars that have lost most of their envelope at the end of the RGB phase. Man y sdB stars 



have a close stellar companion that went through CE phase. This explains the envelope ejection ('e. g..lHan et al 
2003 ). However, many others might be formed through a CE phase with a substellar companion (ISokei]ll998h " 



Although some substellar companions might survive the CE phase, th ere are recent indications that many of 
the substellar companions merge with the RGB core (jCeier et al.ll201l[ ). 



Ultraluminous core collapse SNe. There are hints that so me ultraluminous core collapse SNe exp erience an 
extreme mass loss episode ^ 1 — 100 yr before explosion fe.g.. IChomiuk et al.ll2011t see discussion in IChevalier 
20121 ). I speculate that this coincidence can be caused by a companion that spirals-in inside the envelope and 
collides with the core. The CE process causes the extreme mass loss rate, while the collision with the core 
expedites the evolution toward explosion. The companion can be several solar ma sses main sequence star, or a 
neutron star. A similar idea was put forward simultaneously bv IChevalien (|2012i ). This speculative scenario 
will be studied in a future paper. 
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Complicated planetary nebula structures. The CE merger can l ead to a PN with v ery comphcated structures. 
The PN NGC 5189, for example, has a very complicated structure (jSabin et al.ll2012l ). and future studies should 
examine whether CE merger can indeed explain such a structure. 

Missing long-orbital period post-CE binaries. The migration process at the termination of the CE phase 
reduces the orbital separation even if a merger does not occur. This might explain the finding that all post 
CE binary syst ems in the homogeneous SPSS sample have orbital separation well below expectation of the CE 
a-prescription (jRebassa-Mansergas et al.ll2012l ). 



I thank Amit Kashi and an anonymous referee for helpful comments. This research was supported by the 
Asher Fund for Space Research at the Tcchnion and the Israel Science Foundation. 
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